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I. INTRODUCTION 

The first notion of fractional excitations in condensed 
matter physics goes back to the appearance of soliton 
mid-gap modes in polyacetylenoi, where the effective net 
charge of one kink excitation is e/2, i.e., one half of the 
electron charge. At a similar time, the field of fractional 
statistics, founded in the work by Leinaas and Myrheiro^, 
attained broad attention due to the work by Wilczek in 
19822ii. In strongly-correlated many-body systems, the 
phenomenon of fractionalization, where the elementary 
excitations of the system carry only a fraction of the 
quantum numbers of the constituents, has become known 
to occur in a variety of cases. 

The first physical system in which fractional excita- 
tions and the associated fractional statistics have been 
discussed on a unified footing is the fractional quantum 
Hall effeci^i^i (FQHE). There, the quantum statistics of 
the anyonic quasiparticles can be understood in terms of 
a generalized Berry's phased, which is acquired by the 
wave function as quasiparticles wind around each other. 
This is a sensible concept in two dimensions only where 
one can uniquely define a winding number for the braid- 
ing. In the FQHE, the fractional statistics is known to 
occur in the presence of a magnetic field violating par- 
ity (P) and time-reversal (T) symmetry. In recent years, 
there have been tremendous efforts to study the frac- 
tional excitations of the FQHE experimentally in order 
to confirm the prediction from theory and to validate 
fractional statistics as a concept being realized in nature. 
This, however, has remained inconclusive in certain as- 
pects and thus is still a subject of current discussion and 
work^iiaiiiii^iH. 

Later, the concept of fractional statistics has been 
found to occur in one-dimensional spin-1/2 antiferro- 
magnets, where it can be defined in terms of a gener- 
alized Pauli principle obeyed by the excitationsi^ and, 
as shown recently, by a phase the wave function acquires 
when two spinous move through each other—. The frac- 
tional charge of the quasiparticles in the FQHE corre- 
sponds to the spin 1/2 of the elementary spinon excita- 



tions in these systems, which is fractional as the Hilbert 
space is built up by spin flips which carry spin 1. As 
one-dimensional systems are amenable to a host of exact 
methods, many exactly solvable models exhibiting this 
behavior existji^'^ii^'^. In particular, various properties 
of fractional excitations in spin chains have been observed 
experimentally2£i^i2^. 

In general, it appears to be that P and T violation is 
intimately related to the occurrence of excitations obey- 
ing fractional statistics in two dimensions, which both 
applies to quasiparticles in the FQHE and spinous in a 
quantum antiferromagnet. These symmetries may be ex- 
plicitly broken as in the FQHE or generated by sponta- 
neous symmetry breaking. For two-dimensional antifer- 
romagncts, the concept of fractional excitations is less 
established than for the one-dimensional case^^. In par- 
ticular, finding solvable theoretical models in which the 
phenomenon occurs has been one predominant area of 
research in the field. Significant progress has been ac- 
complished for dimer model a^^'^^ . 

In addition to important questions with regard to 
the general principle underlying fractional statistics, 
two-dimensional spin liquids are of special interest 
with regard to investigation of the hypothesized link 
between fractionalization and high-Tc superconductiv- 
itySSii^. Moreover, in many systems where fractional- 
ization occurs, there is the ambition to use the topolog- 
ical degeneracy contained in these systems for quantum 
computing, where topological information can serve as 
a quantum bit with negligibly small local dccoherence 
rates2^. 

The paradigmatic state for a S ~ 1/2 spin liquid is 
the chiral spin liquid (CSL) introduced by Kalmeyer and 
Laughlin^i^, which is constructed to spontaneously vio- 
late the symmetries P and T, and can be defined on any 
regular lattice including both bipartite and non-bipartite 
lattices. The universality class of chiral spin liquid states, 
and in particular the order parameter and the topological 
degeneracy^, was defined by Wen, Wilczek, and Zee^. 
A CSL state has also been constructed by Yao and Kivel- 
sonS in the Kitaev model^ on a Fisher lattice, i.e., a 
honeycomb lattice of triangles. Recently, a family of non- 



Abelian CSL states^^ has been proposed for general spin 
S", whose wave functions correspond to the bosonic Read- 
Rezayi series of FQH states^. The non-Abehan statistics 
of the spinons has even been conjectured to be a general 
property of spin S antiferromagnets^. 

As in the one-dimensional case mentioned above, the 
spinons in the CSL exhibit quantum- number fractional- 
ization and carry only half the spin of the bosonic spin ex- 
citations in conventional magnetically-ordered systems, 
which carry spin 1. Whereas the spinon appears to 
be the fundamental field describing excitations in two- 
dimensional 5* = 1/2 antiferromagnets in general, an 
effective description by magnon-like excitations proved 
rather adequate for the generic model. The reason for 
this is that the confinement between the spinons is gener- 
ically so strong that the underlying excitation structure is 
mostly suppressed. In the CSL, however, the spinons are 
deconfined. The model is hence ideally suited to study 
fractional quantization of spinons in two-dimensional an- 
tiferromagnets. In spite of its promising properties, for 
nearly two decades since its emergence, the CSL lacked 
a microscopic model where it is realized. 

In this article, we develop an analytical method for the 
construction of parent Hamiltonians for the CSL. The 
method relies explicitly on the singlet property of the 
CSL, as this allows for a spherical tensor decomposition 
of the destruction operator we introduce. From differ- 
ent tensor components, we construct two different parent 
Hamiltonians, which annihilate the CSL and hence have 
it as a zero-energy ground state. One of the Hamilto- 
nians has been presented in a Letter previously^; both 
Hamiltonians contain 6-body interactions. One of the 
key issues we address here is whether the CSL is the only 
ground state of these Hamiltonians. To answer this ques- 
tion, we perform exact diagonalization studies of both 
models for a 16-site square lattice. In particular, we 
introduce an adapted Kernel sweeping method, which 
allows for an efficient numerical implementation of the 
complex and technically cumbersome Hamiltonians we 
investigate. We find that the model we introduced pre- 
viously has indeed the CSL as its (modulo the two-fold 
topological degeneracy) unique ground state. For the 
other Hamiltonian we present, however, we find that the 
CSL is not the unique ground state. Hence only the for- 
mer model is useful for further analysis of e.g. the spinon 
spectrum. 

The paper is organized as follows. In Section [TTl we 
review the chiral spin liquid ground state and its ba- 
sic properties. After outlining the general construction 
scheme for the Hamiltonians in Section IIIIl we formu- 
late a destruction operator for the CSL state in Sec- 
tion IIVI and exploit the spin rotational invariance of the 
CSL state to decompose the destruction operator into its 
spherical tensor components, which annihilate the CSL 
state individually. The proof that the destruction oper- 
ator annihilates the CSL ground state is given in Sec- 
tion |Vl In Section IVIl we introduce a Kernel sweeping 
method to compute the CSL Hamiltonians. We present 
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FIG. 1: The model is defined on a square lattice length L on 
a side such that the total number of sites is given hy N = 
L^. The image shows the lattice for A'' — 16. The shaded 
circles (including the origin) indicate those lattice sites for 
which G{z) = — 1 and the open circles those sites for which 
G{z) = -1-1. The sites on which G{z) = —1 define a sublattice 
with twice the original lattice spacing; the doubled unit-cell 
is shown as the shaded region in the figure surrounding the 
origin. 



the method in detail and emphasize its applicability to ef- 
ficiently compute n-body interactions for finite-size exact 
diagonalization studies. The numerical results obtained 
with this method are discussed in Section IVlIl We con- 
clude this work with a summary in Section [VIIII 



II. CHIRAL SPIN LIQUID 

The CSL was originally conceived by D.H. Lee as a spin 
liquid constructed by condensing the bosonic spin flip 
operators on a two-dimensional lattice into a FQH liquid 
at Landau level filhng factor v = 1/2. The ground state 
wave function for a circular droplet with open boundary 
conditions, on a square lattice with lattice constant of 
length one, is given hy^^^ 
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where M = N/2. The z's in the above expression are 
the complex positions of the up-spins on the lattice: z = 
X + iy, with x and y integer. G{z) = (— l)(^+i)(!^+i) is a 
gauge factor, which ensures that ([1]) is a spin singlet (see 
Fig. [1]). Lattice sites not occupied by z's correspond to 
down-spins. 

For our purposes, it is propitious to choose periodic 
boundary conditions (PBCs) with equal periods Li = 
L2 = L, L even, and with N ^ L"^ sites. Following 
Haldane and Rezayi^, the wave function for the CSL 
then takes the form 
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where = —'di{—w) = '&i{w\e~'^) is the odd Ja- 

cobi theta function^. The zeros for the center-of-mass 
coordinate Z ~ Zj must he in the principal region 
< Re(Zi) < L, < Im(Zi) < L and satisfy Zx + Z^^ 
L + iL] the freedom to choose Zi reflects the topological 
degeneracy and yields two linearly-independent ground 
states for the CSL. These states are spin singlets, are 
invariant under lattice translations, and are strictly pe- 
riodic with regard to the PBCs. 

III. GENERAL METHOD 

In order to construct a parent Hamiltonian for the chi- 
ral spin liquid, one first derives the destruction operators 
for the ground state. In our formulation, the destruc- 
tion operators are constructed from a set of operators Uj 
where j — 1, . . . ,N indexes the lattice sites. The oper- 
ators Wj, to be introduced in Section llVI below, are not 
themselves destruction operators, but have the property 
that, acting on the groimd state, they produce a result 
independent of the site index j: tUi \ip) = coj \tjj). There- 
fore, once the above result is established in Section |Vl it 
follows that the diff'crencc of any two of the operators is a 
destruction operator for the ground state: c?y = a;, — Uj. 

In order to construct a sensible parent Hamiltonian, 
one must minimally demand that it be a translationally- 
invariant scalar operator. In order to put the Hamilto- 
nian in this form, it is shown in Appendix |A] that the 
operators may be written as ujj — 17^ + 15^ where flj 
and Uj are vector and third-rank spherical tensor oper- 
ators respectively and where the superscript indicates 
the component in spherical notation. The operators ^Ij 
and 13 j are given explicitly in terms of spin operators in 
Sections T7T\ and ITTbI 

As is discussed in detail in Section llVi the Wigner- 
Eckhart theorem guarantees that all components of the 
operators D^j = fli — flj as well as X>y = I3i — 13 j are 
destruction operators for the chiral spin liquid ground 
state so long as the reducible tensor operator dij is. One 
can then construct Hamiltonians based on either set of 
operators: 

H = Y,^l■^^J (3) 

for the vector destruction operators or 
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for the rank-3 spherical tensor operators. Either Hamil- 
tonian is a scalar and is translationally invariant, both 
of these properties guaranteed by the construction. Ad- 
ditionally, since the Hamiltonians are positive semi- 
definite, the chiral spin liquid is a ground state of the 
model. It should be noted that these models arc not 
themselves unique as one could include any coefficients 
Jij into the sums of Eqs. [3] and 3] and remove the restric- 
tion that only nearest-neighbor sites are summed over. 
These two models do, however, represent the simplest 
models from each class. 

In Section IVIl a numerical method is developed for 
performing the exact diagonalization of these Hamilto- 
nians that can handle the large number of interactions 
efficiently. This method is used in Section IVIII to show 
that the model given by Eq. [3] has exactly two ground 
states, as expected due to the topological degeneracy of 
the chiral spin liquid on a torus, and that these states 
are precisely the chiral spin liquid ground states given 
in Section |IT] above. Adopting the same procedure, the 
Hamiltonian given in Eq. |4] is shown to have a larger 
ground-state manifold which is not exhausted by the chi- 
ral spin liquid ground states. 



IV. ANNIHILATION OPERATOR FOR THE 
CHIRAL SPIN LIQUID 

The Hamiltonian which stabilizes the chiral spin liquid 
is generated by first finding a set of operators LUi, where 
I is a site index. These operators are not themselves 
destruction operators, but the bond operators uji — loj, 
where i and j are any two distinct sites, will be shown 
to destroy the CSL ground state. The operators may be 
written as ujj = lij^ — uJ where lo'^ = Tj + Vj and 

= li2K,,,S;S^(l + S!) (5) 

The two sets of coefficients Uij and Kijk are defined in 
Section IIV CI below and the prime on the sum indicates 
that one must exclude the coincidences of i and k. 

The operator loJ is related to by a 7r/2 rotation 
about the x-axis that maps Sz and Sy into —Sz and ~Sy. 
This means that the entire operator uoj is given by 
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In writing down Eq. [3 the fact that X^i^^j ^ij = 0' been employed. This will be demonstrated in Scc- 
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tion IIV CI below. While the operators uji are not them- 
selves destruction operators for the CSL ground state, it 
will be shown in Section W\ that da = l 



LOj is a de- 



struction operator for the ground state for any choice of 
i and j. 

The operators ujj are reducible and can be decomposed 
into irreducible tensor operators, in this case of ranks 1 
and 3. From Eq. [7] it is clear that every term except for 
the 5? SI term is the (or z) component of a rank-1 
(vector) operator. This final term can be decomposed 
into rank-3 and vector components. 

It is straightforward to show that if an operator d is a 
destruction operator for the CSL ground state, then each 
of its irreducible components are as well. This is because 
the Wigner-Eckhart theorem tells us that acting with an 
operator T^^ on a state \nqmq) with angular momentum 
q and z-component niq gives 
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where n and n' are any quantum numbers other than 
angular momentum. Since the CSL is a spin singlet: 
q = ruq = 0, it follows that there is only a single non- 
zero term in the above sum corresponding to j' = j and 
to' = m. This means that by decomposing the destruc- 
tion operator for the ground state d into its tensor com- 



ponents, which may be written d = aj Tq, acting on 
the ground state to obtain 



Q^d\^)=Y,aj I"'i0) , 
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and noting that states with different values of j are neces- 
sarily orthogonal, it immediately follows that each of the 
states in the sum are themselves zero and hence the op- 
erators are destruction operators for the ground state. 
In Sections IIV Al and IIVBI we give two classes of opera- 
tors that are obtained from the reducible tensor operator 
Uj in Eq. [71 



A. Vector destruction operator 

As shown in Appendix El the operator Sf Sj may 
be written as the sum of the 0-components of a vector 
and a third-rank tensor. The vector component is given 
by 



- [(5. • Sj) SI + {S, ■ Sk) St + {Sk ■ S,) 5|] (10) 
and. working from Eq.[71 the vector operator flj is given 

by ' 
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Since — flj is a destruction operator for the ground 
state, it immediately follows that one may construct a 
Hamiltonian for which the chiral spin liquid is the exact 
ground state as 



positive semi-definitencss in H) in constructing a parent 
Hamiltonian for the CSL. 



ij = ^ (r2j - Oj)^ • {n, - 



(12) 



B. Tensor destruction operator 



where the sum runs over all nearest-neighbors on the lat- 
tice. By construction, the Hamiltonian is a scalar opera- 
tor and translationally invariant. 

However, note that there is nothing restricting possible 
models to run only over next-nearest neighbors. Rather, 
one can consider any combination of bond-operators (in- 
cluding arbitrary coefficients so long as one maintains 



It is also possible to create a set of third-rank tensor 
destruction operators. As shown in Appendix El the 
operator Sf. may be fully decomposed into the 0- 

components of a vector operator (given in Eq. llOp and a 
third-rank tensor operator, which is necessarily just the 
difference between Sj and the operator in Eq. [TOl 
This gives a destruction operator whose 0-component is 
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The other eomponents are straightforward to obtain (see 
Appendix E| and one may again use these operators to 
form a Hamiltonian for the chiral spin hquid according 
to 
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The Hamiltonian in Eq.[T3]has two significant advantages 
over the model in Eq. 1121 it depends only on one set of 
coefficients {Kijk but not Uij) and, because the operator 
in the sum in Eq. [13] is symmetric under interchange of i 
and fc, one may replace Kijk by A^jk = {Ki-ju + Kkji) /2 
where the new coefficients are manifestly symmetric in 
interchange of the first and third indices. Unfortunately, 
it turns out that the CSL is not the only ground state of 
this model, as will be discussed in detail in Section IVIII 



Coefficients 



The coefficients appearing in Eq. [7] are functions of 
the distance between the sites of the form Kijk = 
K{zk — Zj, Zi — Zj) where 



K{x,y) 



N/2 



^— lim y 



0<zo<-R 



P{x - ZQ,y) 
X - zo 



(15) 



and the sum over zq is a sum over all lattice translations: 
^0 = (w + in) L for m and n integer. This sum guar- 
antees that the function K{x,y) is periodic in its first 
argument. 

The coefficients Uij = ttU{tt [zj — Zi] / L) / L are given 

by 
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where W{z) is the periodic extension of 1/z to the torus^ and also related to the logarithmic derivatives of the theta 
fimctions: 
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The function P(x, y) is given by 

P{x,y)= YwTv y 



Co(^ [zo-y]) e-l^l^«-^l' 



(18) 



where Go{x) = cos a: + cosh x and where n{y) is a normal- 
ization factor chosen such that P(0, y) = 1 which entails 
the choice 



n(z) = z93(^Rc[z] 
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While the form of the coefficients as given by Eqs. [T5] - 
[17] are essential for forming a Hamiltonian that stabilizes 
the CSL, there is significant freedom in how one chooses 
the function P{x,y). The only requirements arc that 
it be a periodic function of y, fall off faster than 1/x 
with increasing a;, and be analytic apart from first-order 
poles that occur at the coincidence of the two arguments: 
X = y. It is straightforward to show that U{z) is an odd 
function; this in turn guarantees that Uij = and 
lets this sum be dropped, as was done in writing down 
Eq.[7] 



V. PROOF OF SOLUTION 

In order to prove that cither of the Hamiltonians given 
in Eqs.[3]and[4]are true parent Hamiltonians for the chiral 
spin liquid, we must demonstrate that ujj \ip) = uji j?/;) 
which we will demonstrate by first showing that 



{zi-- ■ ZM\uJj \xp) = /(Z) (zi • • • ZM\i>) 



(20) 



where f{Z) is a function only of the center of mass: Z = 
2i- This identity in turn follows from the fact that 
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(21) 



and the result that the function f{Z) is both odd and 
periodic. To see this, recall that one can write ujJ = 
U uj'j U'^ where U performs the 7r/2 rotation about the 
X-axis as discussed in Section lTVl abovc. The CSL ground 
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state is invariant under such a rotation so that 

(zi •••ZM |c^+|^> = (zi---ZM|(7Wi7|V'> 
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where {wi}, the locations of the down spins on the lattice, 
is the complement of {zi}. It then follows from Eg. [2T] 
that 
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Assuming that the origin of the lattice is chosen such 
that the sites occupy positions Zi ^{i + im) for i and m 
integer, it is straightforward to show that 



L(L-l) , 
Z + W^ \ ' {l + i) L, 



(24) 



and since L is even it follows that the sum of Z and W is 
equivalent to a translation of the lattice zq. Because the 
function f{Z) is periodic and odd, both properties will 
be shown below, it immediately follows that /(W) = 
/(zq — Z) = —f{Z). Combining this fact with Eq. [23l 
completes the proof that Eq. [21] entails Eq. [20l 



Action of Ti 



In order to prove Eq. [21] we first consider the off- 
diagonal terms in the operator oj^ which come from Tj 
defined in Eq. [3 We consider a general element of the 
vector Ti lib): 



(zi---zm|Tj IV') = X XI -^^J*^ 
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The element is clearly zero unless Zj e {zi • • • zm}- When this is satisfied, acting onto the bra on the right-hand side 
of the equation with the spin operators wipes out the matrix element unless z^ S {zi • • • zjv/} and replaces Zj with z^: 
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The upper limit of M = N/2 (rather than N) on the sum on i indicates that z^ must be a member of the up-spins. 
Rewriting Kijk — K{zk — Zj, Zi — Zj) and defining z ~ Zk — Zj, this may be rewritten as 
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Using the definition of the coefhcicnt K from Eq. [151 this can be rewritten as 
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Since the wave function itself is periodic, the two sums over z and over zq may be combined into a single sum that 
runs over the entire infinite lattice for which we use the variable x = z — zq. However, since the point z = is missing 
from the original sum, all of its images in the infinite lattice will be missing from the second sum and this must be 
subtracted off, giving 
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Dividing both sides of the equation by (zi • • -zmIV) ^-nd rewriting the ratio of elements in terms of the analytic 
function of x, A{x) given in Appendix [Cl vields 
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Note that A{x) is an analytic function only of x, and not of the remaining {zi} on which it also depends. 

The first sum in Eq. [30] may be evaluated with the corollary to the Singlet Sum Rule. Eq. IB13I A derivation of 
the sum rule and the necessary corollary is given in Appendix [BJ The function P{x,y) falls off exponentially with 
increasing x while the quantity A{x) G{x) is essentially constant due to the periodicity of the wave function. 

This guarantees that the sum is absolutely convergent and the sum rule may be applied. Additionally, the product 
P{x, Zi — Zj) A{x) is itself an analytic function of x. As a function of x, the function P{x, Zi — Zj) necessarily has 



poles. However, these occur when x = Zi 



and the function A{x) has second-order zeroes at these locations since 



this corresponds to a coincidence of up-spins. Since the product is analytic and the sum is absolutely convergent, the 
singlet sum rule may be applied to give 
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Using the fact that A{Q) = P{0, Zi — Zj) = 1 and the relation for dA/dx given in Eq. IC9[ this becomes 
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The sum on i may be completed for the terms containing the W functions (picking up a factor of M = N /2 — 1) and 
this gives, renaming ^ as i. 
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where 



y—1 



(35) 



The fact that f{Z) is both odd and periodic, required 
for the proof of Eq. [^D] above, follows from these same 
properties of the W function. Comparison with Eq. 1161 
shows that 



(zi ■ ■■zm\Tj 
(zi • • •za/IV') 



M 



/(Z)-^f/,, 



(36) 



if Zj is an element of the up-spins and zero otherwise. 



The matrix element vanishes unless both Zi and Zj are 
elements of {zi • • • zj\/}. Therefore, the diagonal contri- 
bution to the operator ujj gives 
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if Z; 



G {zi • • • za/} and otherwise. Combining Eqs. 
and [36] yields Eq. [21] and therefore proves that the chi- 
ral spin liquid is an exact ground state of either of the 
Hamiltonians in Eqs. [T^ or [HI 



VI. KERNEL SWEEPING METHOD 



B. Action of Vj 

The action of the operator Vj on the CSL ground state 
is straightforward to compute. Proceeding in an analo- 
gous manner, we have 

(zi • • • ZmI Vj Wl = 

E • • • ^^^1 + ) + ^i) I'/') • (37) 



To implement the Hamiltonians given in Eq. [T^] and 
Eq. [M] one has to take into account that 6-body terms 
appear in the Hamiltonians. For microscopic models con- 
taining many-body interactions, one must be very effi- 
cient if one hopes to write down the Hamiltonian in a 
reasonable amount of time. For our Hamiltonians, this 
is because there are, even for a lattice with only = 16 
sites, literally thousands of terms in the Hamiltonian cor- 
responding to all the different ways to choose six sites out 
of sixteen. In contrast, a model with only two-site inter- 
actions on the same lattice would only have 15 terms to 
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compute after taking into account translational symme- 
try, even if the model had infinite range. In this section, 
we describe an algorithm for calculating the Hamiltonian 
very efficiently, called the kernel sweeping method. 

As an example to illustrate the kernel sweep- 
ing method, we will consider the computation of a 
Heisenberg-type Hamiltonian such as 



H — ^ ^ Jij Si ■ Sj 



(39) 



where, since we are dealing with a Hermitian operator, we 
only need to include the upper triangle. The extension 
of this array to a p-site operator is straightforward; in 
that case one must consider the action of the operator 
on a 2P-dimensional basis. Therefore, the corresponding 
operator for the chiral spin liquid Hamiltonian given in 
Equation [12] is 64-dimcnsional. 

The code next loops over all possible values of i and j 
and docs the following. First it computes 



We work in an basis and label the states by a binary 
number where up-spins are treated as 1 's and down-spins 
are treated as O's. We first note that since this is a two- 
site interaction, in order to implement this model all we 
really need to know is how the operator Si ■ Sj acts on 
the four-dimensional basis \siSj). This action may be 
summarized as 



Ri 



{2-1, 2^-1}, M„.{2-\ 2^-1}] 
„„} . (43) 



All this means is to compute the contribution of the two 
spins at sites i and j to the binary number that will label 
the entire state. For our example, assuming that we are 
at a point in the loop where z = 3 and j = 7, this gives 



(iil 


ati 


(Ul 


(TTI 




1/4 











lii) 





-1/4 


1/2 





lit) 





1/2 


-1/4 





in) 











1/4 


ITT) 



(40) 



where the table format shows the order of the basis vec- 
tors. It is only necessary to compute this matrix once 
at the beginning of running the code. One stores this 
matrix as a set of rules 



(41) 



where the {s}„j and the {s},^ are a binary shorthand for 
the states in this two-dimensional basis and ^Imn sltc the 
elements in the matrix. In this example we would have 



R = <i [00,00] ^ i, [01,01] 



[10, 10] ^--,[11, 11] 



-i,[01,10] 



(42) 



i?37 = {[0,0] ^:^,[2^26]^-:|I, [2^,2^]^:^, 

[22,22] [22 + 2^,22 + 26] ^:|l|.(44) 



The code next computes the contributions to the binary 
numbers labeling the states from all the sites that are not 
involved in the interaction. There are 2^~p of these and 
for our two-site example this list is 



N 



Finally, one updates the Hamiltonian according to 



(45) 



H = H + Rij® Bij 



(46) 



where the addition means to add the matrix defined by 
these rules and the generalized outer product means 



Ri 



'B^J = {[{s}^ ■ {2-1,2^-1} +6, M„ • {2-1,2^-1} +&] ^ J., a„„} 



(47) 



for b an element of Bij . In this way one may construct the 
Hamiltonian extremely quickly since all the steps involve 
list operations and there is only a single loop over the 
N choose 2 ways to pick the sites i and j. (In practice, 
one uses translational invariance to fix i = 1 and, for a 
two-site operator as in this example, the loop is then over 
the N — 1 ways to choose the remaining site.) 

Let us now work this out explicitly for the vector 
Hamilton operator Eq. [T^l Setting Di, — fli ~ , we 



split up the Hamiltonian into 



(ij) 



(48) 



where the z-component as well as the ladder components 
of the vector operators can be written out in terms of spin 
operators S'^ 5+ = 5^ -f iSy, and 5" = 5*^ - iSy. As 
the treatment is very similar, we constrain our attention 
to the contribution J2{ij) ^tj'^^tj^ where for clarity we 
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again write out the + ladder operator explicitly: 



E 



K 



ijk 



i^j^k - Sj'Sk) + - {Sj ■ Sk) 5;* 



--(Sk 



Si) Sj' 5 ' 



St 



(49) 



Using the notation analogous to Eq. [U 

TV 1 f TV 



B^^{Y,sa'-'), (50) 



we can write 



E 



E ^^Ok 
k 



B}jk 



+ Rf ® Bf 



(51) 



where R^-^ and i?" relate to the first and second sum 
of Eq. I49i respectively. Given these 3-body operators in 
above notation, the total 6-body interaction can be con- 
veniently computed. The implementation of the tensor 
Hamilton operator Eq. [T3]is completely analogous. 



VII. NUMERICAL CONFIRMATION 




X 



K 



M 



FIG. 3: Low energy spectrum of the Hamiltonian Eq. 1121 
scaled down to order of unity. There are two E = eigenval- 
ues at the r point. 



Using the method outlined in Section I VII above, the 
models in Eq. [T2] and Eq. [T3] have been solved by exact 
diagonalization on 16-site lattices with periodic boundary 
conditions. We start by considering the vector Hamilto- 
nian given by Eq. 1121 The spectrum is shown in Fig. [3l 
the points in the Brillouin zone which label the axis of 
this figure are shown in Fig. [2l We find the spectrum to 
be positive semi-definite, with a doubly-degenerate zero- 
energy state at the F point. The rest of the spectrum 
is well separated from the ground state by a gap that is 




•X 



FIG. 2: A plot of the symmetry points in the first Brillouin 
zone. The arrows show the path taken in plotting the energy 
spectra in Fig. |31 starting from the origin at F = (0, 0). 



substantial and we believe not due to finite size effects in 
the calculation. This claim is based on the fact that it 
exceeds the finite size level splitting of the spectrum by a 
factor of ^ 15. The presence of a gap is expected between 
the chiral spin liquid ground state and what should be a 
two-spinon excited state. The spinon excitations of this 
model will be addressed in future work. 

We now discuss the two orthogonal zero-energy eigcn- 
states. For comparison, we construct the CSL state Eq.[2] 
explicitly and find a two-dimensional subspacc of func- 
tions with the center of mass variable being treated as 
an external parameter. We have computed the overlap 
of the Hamiltonian ground state subspace and the CSL 
subspace and find that they match perfectly. Therefore, 
the ground state of this Hamiltonian is indeed the two- 
fold degenerate CSL state. Additionally, we have only 
two zero-energy states, by which follows that the CSL 
state is the only ground state of the model, a statement 
which cannot be achieved analytically. 

For the tensor Hamiltonian, however, we find that 
the zero-energy subspace is massively degenerate. It of 
course contains the CSL, in accord with the analytical 
proof, but also many additional states. While the restric- 
tion to small system sizes prevents us from studying the 
thermodynamic limit precisely, our numerical findings in- 
dicate that the Hamiltonian Eq. [T3]does not stabilize the 
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CSL state as the unique ground state, which thus singles 
out the model in Eq. [T2]to be subject of further study. 

VIII. CONCLUSION 

In this work we have shown a method for constructing 
parent Hamiltonians for the chiral spin liquid. We have 
computed the spectra of the Hamiltonians by use of a 
Kernel-sweeping method in exact diagonalization. There, 
for the Hamiltonian operator composed of the spherical 
vector component of the CSL destruction operator, we 
observe that the CSL states are the only ground states 
of the model. We conclude that this model is a promising 
candidate to also study the elementary excitations of the 
model, i.e., spinous, and many other questions in the field 
of two-dimensional fractionalization of quantum numbers 
in spin systems. 
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APPENDIX A: TENSOR DECOMPOSITION 

The operators ojj introduced in Section ITVl mav be de- 
composed into irreducible spherical tensors of ranks 1 



and 3. We write these irreducible operators as T^; q and 
m correspond to angular momentum and its z compo- 
nent respectively. We wish to write to = "^CqT"^, where 
is the collection of all operators which transform as a 
spherical tensor of rank q. Here we have suppressed the 
site index on the operator lo. 

The operator in Eq. [7] that is not manifestly the com- 
ponent of a vector is Sj S^, which is a component of 
a third-rank Cartesian tensor. In order to keep the nota- 
tion manageable, we start by considering the direct prod- 
uct of two operators U and V with angular momentum 
ji and j2 respectively. An element in the direct product 
space of these operators may be written as 



u?,i,v^^= E c'-^t^- (A1) 

jl2 = \j2—jl\ ™12 = — il2 



in terms of irreducible spherical tensors T^^^ carrying 
angular momentum ji2 with z-component toi2 — mi + 
1712. Eg. lAll mav be inverted to give 



nii=-ji m2=-j-2 

Using these equations, one may construct correspond- 
ing expressions for the product of three vector operators 
by applying Eg. lAll twice: 



TJ31 yj2 T^J3 

^ mi 7712 m3 



ii+j2 ii2 

Z-^ il i2 312 mi2 '^'nii 

jl2 = -\jl-j2\ mi2 = -jl2 

jl+j2 jl2 jl2+j3 j 

E\ ^ ^mim2mi2 \ ^ \ ^ /^mi2m3m /Tij(ji2 

il i2 jl2 312 33 j in 

3l2 = -\3l~32 \ mi2 = -jl2 j=\ji2-j3 \ m=-j 



(A3) 



The second superscript on the tensor T in the last line distinguishes between the different tensors of the same rank 
that appear when combining three vector operators; since 1®1(3)1 = 3®2©2©1©1®1®0, there are two rank-2 
spherical tensors and three vector operators that can be formed. For the case of interest mi = m2 = — and 

il = j2 = js = 1 , the expression reduces to 

U^V^W^= E E' C^i2Cl2^n 

jl2=0 j = |ji2-l 

which shows that the operator contains only vector and rank-3 tensor components, but no scalar or rank-2 tensor 
components. Note that the second index on the rank-3 tensor has been suppressed since the construction of this 
object is unambiguous. 



1 rpHO) 

'V3 ° 



1(2) 



' 



(A4) 
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Applying Eg. IA2I twice, the rank-3 tensor component is 

2 1 

mi2= — 2 m3 = — 1 
1 

'-'2 1 S'-'I 1 2 "^imi l^ma '^^^ma l^Oj 

mi ,7712 5^3 — ^1 

5£/^ ~{V - W) U"" -{W -U) V' ^^^^ 



10 

where we have used the fact that the dot product is U • V = X]jti(^1)™ ^-m i^i the spherical representation. A 
similar construction can be used to find the vector operator or, one may note from Eqs. IA4I and IA7I that the vector 
component is equivalent to 

1^- - ■ ^) + ■ ^) + ■ ^) (A8) 

as used in writing down Eq. 1111 

Construction of the remaining (a:: and y) components of the vector operator in Eq. IA8I is straightforward since one 
merely replaces z with either x or y. In order to construct the remaining six components of the rank-3 tensor operator 
one simply applies Eq. IA2l twice without specifying m = 0: 

1 

rp3 _ \ ^ ^771-77737713771 ,^rni7ra2"7-ni3 7-7-1 T/1 T/t/1 / A Ct^ 

L"2 1 3 "-"l 1 2 "^7771 '^777 2 ^''^7773 ' > 

7771 ,7772,7773 — — 1 

The explicit form of these components are 

T'l = ^ \hV''W'' ~V -W) U+ +{5U'W^ -U -W) V+ +{5U'V' - U -V) W+] (AlO) 

2 v30 



T| = [U+V+W'- + U+V''W+ + U''V+W+] (All) 

2 v3 

r| = --l^u+V+W+, (A12) 

with the remaining three components obtained from T^^^ = (—1)™ (T^^J^. 

I 



APPENDIX B: SUM RULE 

The sum rule used in Section |Vl on which the proof 
that Lo destroys the ground state hinges, is given by 



hm y 

Q<\z\<R 



G{2 



0. 



(Bl) 



The sum rule has been first stated in a mathematical 
framework of coherent state systems by Perelomov and 
has later been re-derived by Laughlin^i^; in this ap- 
pendix, we show how to obtain the sum rule by appli- 
cation of Jacobi's imaginary transformation. We first 
consider the related sum 



F(c) 



lim 

J?-»0 



G{z) exp 



0<|2|<fl 



1 



TT I ,2 

2 



(B2) 



In order to prove the sum rule in Eq. IBU we will first 
show that -F(c) = for any value of the parameter c and 



then use this to prove Eg. lBll bv taking derivatives of the 
function F{c). 

In order to show that F{c) — 0, we use the gauge 
function G{z) to write Eg. IB2l as two sums, one over the 
entire lattice and one over the points z' on the lattice 
for which G{z') = —1. As shown in Figure[2l these sites 
defines a sublattice with twice the original lattice spacing. 



F{c) 



62 ' 



2 ^e5=^ 



(B3) 



Setting z' = 2 z we can write this as 
F[c) =^65==^-^ 1^ 



2 



(B4) 



where both sums now run over the entire lattice. Writing 
z = X -'r iy this function can be factored into four sums 
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over the integers x and y: 



F{c) = 



,i(c2;-7rx^) 




y 



In terms of the third Jacobi theta function^S 



E 



(B5) 



(B6) 



this function may be recast as 



F{c) 



2i 



(B7) 



The fact that the two terms in this expression precisely 
cancel is a result of Jacobi's imaginary transformation^, 



(B8) 



and the fact that the third Jacobi theta function is even. 
Application of this identity to either product of theta 
functions in Eq. IB7I shows that the two terms precisely 
cancel, proving that F{c) = 0. This in turn proves the 
n = case of Eg. IBll bv simply setting c = 0. The other 
instances of the sum rule are obtained by noting that 



1 d" 
m\ dc' 



:nc) 



lim > 

R^oo ^ 

Q<\z\<R 



G{z) z"e""e-5 1^1 



(B9) 



where A{z) is an analytic function of z. Since it is ana- 
lytic, we can expand the function A(z) in a Taylor series: 



1^ 



(Bll) 



and, so long as the sum in Eq. IBIOI is absolutely con- 
vergent, we can interchange the order of the two infinite 
sums to obtain 



V i — 

^J\~d? 



lim > ; 

0<|z|<_R 



-^G{z) e 



(B12) 



All terms for which i > 1 immediately vanish from the 
interior sum due to the sum rule in Eq. IBll The term 
with £ = also vanishes because in that case the sum- 
mand is an odd function summed over the entire lattice. 
Finally the term with £ = 1 can be evaluated using the 
sum rule and is simply the negative of the value of the 
summand at z = (which is not included in this sum 
but is included in Eq. IB1|) . Therefore, so long as A{z) is 
chosen so that the sum itself is absolutely convergent. 



lim 



0<|z|<fl 



G[z) 



dA 

dz 



(B13) 



Since F{c) is for all values of c, setting c = in the 
above expression gives the desired result in Eq. IBll It 
should be noted that this proof can be generalized for 
arbitrary lattices. However, for lattice structures with 
non-orthogonal vectors spanning the unit cell, a decom- 
position into Jacobi Theta functions is not possible any- 
more and to follow the above line of proof one has to ap- 
ply the generalized Liouville theorem for two-dimensional 
Riemann theta functions^. 



a. Corollary to Sum Rule 



APPENDIX C: THE FUNCTION A{z) 



The ratio of wave function coefficients appearing in 

Eq.m 



{zi---Zj +x---zm\ 



■ ■ Zm\ 



(CI) 



We now consider the case where we wish to evaluate a 
sum of the form 



lim - A{z) G{z) e 



R-^oa '■ — ' Z 
0<|z|<_R 



(BIO) 



can be written in terms of the Gauge function G{x), the 
Gaussian 1^' , and an analytic function of a;, A{x). 
To see this we note that this ratio may be written explic- 
itly as 
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JJ^ ^(f [Z - Z,]) y ^^{i [z, - z,]) G(z,) G{x) ,f ,f [.-I.P] ^ ^ 

This simplifies by noting that the exponential terms obey an addition formula 



(C2) 



■ [{Zj+xf-\Zj+x\'^] 



e2 



and the Gauge function obeys an addition formula given by 

G{z,) G{x) 



(C3) 



(C4) 



Since the terms involving x* cancel on multiplying the two expressions in Eqs. IC3 land lC4| the ratio of coefficients is 



{zi - ■ ■ Zj + X - ■ ■ zmH') 
{zi - ■ ■ Zj ■ ■ ■ zm\iP) 

where A{x) is an analytic function of x: 



= -A{x) G{x) e 



The derivative of this function is given by 

Evaluating this at x = and noting that Aifi) = 1 from Eg. IC6I gives 

2 , J 



(C5) 



(C6) 



(C7) 



dx 



iy=l i^^j ""'J 

In terms of the function W{z) introduced in Eg. 1171 this may be written as 

2 M 



AT- 



L 2L 



(C8) 



dA 

dx 



v=l 



(C9) 



The final expression follows from the fact that the center of mass zeroes are constrained to satisfy 'Y^^ = a,a 
pointed out in Section [Til 
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